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CN ■ We present a mechanism in which models with two Higgs fields can undergo a two stage 

G\ ■ 

phase transition as the temperature falls. The first stage is a conventional second order (or 
weakly first order) transition in which the symmetry is broken. Shortly thereafter follows 

P 



a first order transition with barrier penetration, bubble production and loss of thermal 
equilibrium. For Higgs potentials with CP violation, we show that the second stage of the 
transition has all the features required for weak scale baryogenesis. 
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1. Introduction 

The SU(2) L x U(1)y symmetry of the standard model is broken in the early universe as 
the temperature falls below the weak scale. The order parameter for the phase transition 
is the real component of the neutral Higgs field, which acquires a non-zero vacuum 
expectation value below the transition temperature. The precise details of the transition 
are still the subject of controversy, though the consensus seems to be that it is weakly first 
order for Higgs masses greater than the experimental lower bound That is, at the 

transition, the Higgs vacuum expectation value (VEV) makes a jump from zero to some 
very small fraction of its eventual zero temperature value. 

Two Higgs models have also been investigated recently |3],j5| , though the whole range 
of parameter space has not been fully explored. Previous work has concentrated on first 
finding the zero temperature absolute minimum of the effective potential and then, in much 
the same way as has been done for the single Higgs model, studying the evolution of this 
one minimum from high to low temperatures. For certain combinations of parameters, one 
loop effective potential calculations suggest that the transition can be first order || 

In this paper, we make use of the large number of parameters in a two Higgs model, 
and construct a potential in which the absolute minimum just after the phase transition 
bears no relation to the absolute minimum at zero temperature. At the transition point, 
the VEV is lured the "wrong" way, into a minimum which is only an absolute minimum for 
a short time close to the critical temperature. As the temperature continues to fall, it soon 
becomes a false vacuum, because a second minimum begins to appear whose depth quickly 
overtakes that of the first. At this point a tunneling process will occur with production of 
bubble walls and loss of thermal equilibrium, thus producing the main characteristics of a 
first order transition. As we shall see, it is not difficult to ensure that the discontinuity in 
the VEV is rather large. In fig. |l|, these various stages of the transition are indicated. 

The electroweak phase transition, strictly speaking, is second order (or at best weakly 
first order), and the first order phase transition occurs somewhat later. However, the 
jump in the VEV during the second stage can be large and thus this mechanism could be 
useful in generating baryons at the electroweak scale according to the scenarios recently 
proposed by Cohen, Kaplan and Nelson in They suggest that top quarks, reflected 

from a bubble wall in a CP violating manner, produce a hypercharge asymmetry which is 
converted into a baryon asymmetry in regions of false vacuum where baryon violation is 
rapid. This asymmetry is swept into the true vacuum when the bubble wall passes and, 
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in order that it be preserved, sphaleron processes in this region must be highly suppressed 
by a large scalar VEV. 

In section 2, we present the simplest example of the two stage transition in a theory 
with just two scalar fields. In section 3, we apply the same ideas to the more complicated 
example of the standard model with two Higgs doublets. In section 4, we show that CP 
violation can be introduced into the potential whilst retaining the other important features 
and thus that it is suitable for weak scale baryogenesis. 

2. The Simplest Example 

Let us consider a theory of two real scalar fields fa and fa with two discrete symmetries 
fa — > —(f) i and fa - * —fa- We will think of fa as the standard Higgs, and fa as an 
additional Higgs. The Lagrangian is constrained to the form: 

£ = ^fad^fa + ^fad^fa - V(fa, fa) , (2.1) 

where 

V(fa, (f> 2 ) = k 1( f>t + Ml + k 3 <f>i<f> 2 2 - 2/x?0? - 2 M |0| • (2.2) 

At zero temperature we will arrange the VEV to lie purely in the (f>i direction and to have 
a given value, in which case (f>\ = n\/k\ = (246 GeV) 2 for instance, and the value of the 
potential here is V\ = —[x\jk\. We can also arrange for there to be a second minimum, 
purely in the fa direction at fa = [i^jki with V% = — ^2/^2, which, if it is not to be the 
absolute minimum, must satisfy | V2 1 < | V\ \ . From the second derivatives of the potential 
at the absolute minimum it is easy to read off the zero temperature mass spectrum : 
m\ = 8(if; m| = 2{k^/kx)^\ — 4n?>- Anticipating the generalization to the real world, we 
suppose that there are some "experimental" constraints on these masses, let's say, that 
mi, m 2 > 50 GeV. 

Next, we wish to make sure that the minimum which appears first during the transition 
is the "wrong" one, i.e., the one in the fa direction. To do this we must recompute the 
positions of the minima, this time as a function of temperature. To see the effect that 
we have in mind, it will not be necessary to go beyond leading order in high temperature 
corrections and for this reason we use a simple method, put forward by Linde J[J, to 
calculate the locations of the minima. We compute thermally averaged equations of motion 
and assume that the particle masses are sufficiently light in comparison to the temperature 
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T that we can use the result (4> 2 ) = T 2 /12 as a good estimate of the fluctuations in the 
scalar fields. In this way we obtain the following pair of equations: 



dV 

Wi 
dV 

d(p 2 



= 4k l( pl + 0i [aiT 2 + 2k 3 (P 2 2 - 4^1} ; 

(2.3) 

= Ak 2 <Pl + 4>2 [a 2 T 2 + 2ks4>l ~ 4/il] ; 



where a.\ = k\ + ks/6 and a 2 = k 2 + ks/6. The non-zero (pi and <p 2 minima appear when 
the temperature drops below 

T C1 2 = M and T C2 2 = M, (2.4) 
ai a 2 

respectively. In order that the VEV occur first in the (p 2 direction, we must make sure that 
T Cl < T C2 , even though we will end up arranging them to be not too different. Equations 
( |2.3|) can have a third solution in which both (pi and <p 2 are non-zero. In the case of interest 
to us, however, this will be a saddle point rather than a minimum of the effective potential. 

Below these critical temperatures in this first stage of the transition, then, the minima 
in the effective potential lie at: 



2 _ 4^f - aiT z 
l — 



Aki 

and 



; (2.5) 



<Pi = , <p 2 = ^7 • (2.6) 

As they develop, at some temperature T t , the depths of the two minima become equal and, 
below this temperature, ( |2.5| ) rather than ( |2.6| ) gives the location of the true vacuum. We 
assume here that bubble nucleation and conversion to the low-temperature vacuum region 
occurs shortly thereafter. T t , thus, signals the onset of false vacuum decay and the second 
stage of the transition. The tunneling temperature can be calculated by substituting 



equations (|2.5| ) and (|2.6j ) into the temperature corrected version of ( pT2|) which is: 

Vfa, 02, T) = Vfa, <h) + T 2 /2 [aicPl + a 2( p 2 2 ] . (2.7) 

The values of the potential at the two minima are always given by Vi = —ki<p\ and 
V 2 = —k 2 cp\ respectively, and so, if we want a large jump in the size of the VEV during 
the jump, we should pick the ratio k 2 /k\ to be large, since, when Vi = V 2 , the ratio of the 
VEV's goes as the fourth root of this quantity. This ratio is restricted, however, by the 
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requirement that the mass of the physical Higgs (coming from ki) be not too small, and 
the desire to avoid extremely non-perturbative couplings. The tunneling temperature is 
given by: 



It is not hard to find fc's and /u's which satisfy the three requirements of large /C2//C1, 
T Cl < T C2 , and \V±\ > | V2 1 at zero temperature. For instance, suppose we choose ki = 0.01, 
then 111 = 25 GeV to give a Higgs mass of m\ = 70 GeV. To maximize &2/&1, we choose &2 
at the edge of the perturbative regime, say ki = 0.8. A large value for /C3, such as ks = 0.6, 
will help us satisfy the temperature constraint. The value of \i2 is then constrained by 
the temperature and true minimum constraints; a value of \i<i = 72 GeV will suffice. The 
critical temperatures work out to T Cl = 151 GeV and T C2 = 152 GeV, and the mass of the 
second scalar is m 2 = 233 GeV. The jump takes place at Tt = 139 GeV, at which point 
the VEV changes discontinuously from 32 GeV in the <p2 direction to 97 GeV in the <pi 
direction. 

3. The Two Higgs Doublet Standard Model 

We would now like to consider a more realistic possibility, where we include Higgs 
doublets $1, $2 instead of singlets, and we will include fermion couplings as well. The 
symmetries 0j — > — (pi will be promoted to <&i — > — and we will allow such terms to be 
softly broken. We furthermore assume that only one of these two doublets (specifically $1) 
is responsible for fermion masses, or at least, for the top quark mass; such an assumption 
is usually necessary to avoid flavor changing neutral currents. The potential with these 
constraints is conventionally written 0-0]: 



T? 



t — 



/,l/2 ,1/2 s. 

(k 2 ai — k x a2) 



(2.8) 



Ai($l$i - vf) 2 + A 2 ($S$ 2 - v 2 2 f 

+ A3[($ t 1 $l-^) + (*2*2-^ 2 2 )] 2 

+ A 4 [($i$l)($ 2 $2) - ($l$ 2 )($ 2 $l)] 

+ A 5 (Re($i$2) - V\V 2 cos£) 
+ A 6 (Im($j$ 2 ) - v x v 2 sin£) 2 • 



(3.1) 
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In this form, one can easily be blind toward the existence of minima other than $j = v\ ; 
$2 = v 2 e 1 ^ (because some of the coefficients above may be negative) and, for this reason, 
we find it more convenient to rewrite the potential as: 



The presence of CP violation is signalled by the presence of n\ and the phase of k^. We 
can redefine the phase of one of the scalar fields to eliminate, for example, the phase of 
fcs; we will assume from now on that k$ is real and positive (a similar phase choice was 
used to eliminate one term in (|3.1|) ). The factor of 4 is present so that when we make the 
substitution: 



the potential closely resembles that of the previous section. 

The main difference from the model of the previous section is the much larger number 
of degrees of freedom to be thermally averaged. We can thus expect the coefficients CKj to 
be larger than before. Further, the Higgs fields are not only self coupled but also coupled 
to the gauge bosons and the fermions, though, of the latter, only the top quark has an 
appreciable Yukawa coupling X t , and this, only to $i- 

We will adjust the phases of the Higgs doublets so that the only fields in ( |3.3|) which 
get a VEV are 0i, (p 2 and ip 2 . Proceeding in the same fashion as before, we can compute 
three thermally averaged equations of motion: 



V($ u $ 2 ) 



4{fci($J$i) 2 + fc 2 (^$ 2 ) 2 + fc 3 ($|$ 2 ) 

+ fc 4 [($l$l)($|$ 2 ) " ($!$ 2 )($2$l)] 

- - ^2^2 - ^Re($l$ 2 ) - /4M$i$ 2 )} • 



(3.2) 





(3.3) 



< 



d(p 2 



dV 



> 







4Mi + [« 2 T 2 + 2(k 3 - k^l + Ak 2 ^ 2 2 - 4^ 2 ] - 2/*|0i ; 








Ak 2 ^l + 1P2 [a 2 T 2 + 2(k 3 + k 5 )(j)j + 4k 2 (/>l - 4^ 2 ] - 2 f 4<f> 1 ; 



(3.4) 



< 



dV 



) 







Ak 1( f>l + </>! [a x T 2 + 2{k 3 - k 5 )<p 2 2 + 2{k 3 + k^\ - 4//?] 
- 2/x§0 2 - 2\i\i\) 2 ; 
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where 

2fcl + ifc3 + i fc4 + ^( 1 + f° s2 M +2A g. 



sin ^yvK /g 5^ 

2, i, e 2 (l + 2cos 2 M 
«a = 2^2 + |fc 3 + 1*4 + \ 2 , 

sin zt7VK 

where e and are the electromagnetic coupling and electroweak mixing angle respectively. 
For the moment, let us restrict our attention to the case when k± = k$ = and ^3 = ^4 = 
so that there is no CP violation in the potential and the VEVs may both be taken to be 
real. With this simplification the equations ( |3.4| ) become the same as ( |2.3| ) except that, 
as expected, there are more contributions to a^. We will take X t = 0.5 corresponding to a 
zero temperature top mass of 123 GeV. 

The analysis of the transition is now almost identical to that described in the previous 
section. The eventual VEV is still <pi = 246 GeV, so if we choose as before k\ = 0.01, 
Hi = 25 GeV, then again mi = 70 GeV. Once more, choosing hi = 0.8, and this time 
&3 = 0.2, we then immediately obtain from (|3.5| ) a>i = 0.987 and a.2 = 2.07. From ( |2.4| ) 
we then obtain T Cl = 50 GeV. The choice \i<i = 43 GeV will then satisfy the remaining 
constraints. We find m 2 = 130 GeV to be the mass of the other neutral scalar as well 
as the charged scalars and the pseudoscalar. They are all equal in mass because of the 
residual symmetry in the Higgs sector when no component of $2 has a VEV. From ( |2. 1| ) 
we find that the first transition occurs at T C2 = 60 GeV and from ( |2.8|) we see that this is 
followed by a tunneling transition at Tt = 46 GeV, at which point, the VEV changes from 
30 GeV in the </>2 direction to 90 GeV in the <f>i direction, as we can see by substitution in 
(|2.5|) and ( p.6| ). Some of the details of this transition are displayed in fig. []J 

Reintroduction of certain of the missing terms causes no substantial changes. The 
main effect of a positive is to increase the charged scalar masses; these do not interest 
us at the moment, so we will leave this parameter vanishing. If we increase h$ to a non-zero 
value, the equations (|321) change, but not their solution. The extra scalar and pseudo- 
scalar masses will also be split, but these are so large that we can allow fcs to be almost as 
great as ks without running into any difficulties with experiment. The main effect of k$ is 
to assure that the phases of $1 and $ 2 are aligned (that is, 1^2 = 0) when they are both 
non-zero, so, even with these complications, we expect the phases of the two Higgs fields 
to be aligned in the transition region comprising the middle of the bubble wall. This term, 
however, will have very little effect on the regions near either of the two local minima, 
when one of the fields vanishes, and it is these regions that constitute the top and bottom 
of the wall. The effects of the (13 an d n1 terms will be studied in the next section. 
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4. CP Violation and Weak Scale Baryogenesis 

We wish, now, to show that the transition which we have described may be useful in 
generating a baryon excess at the weak scale. The potential already has the necessary out 
of equilibrium phase transition, and electroweak interactions can violate baryon number, 
but we have not introduced CP violation. To achieve this, we will now reinstate ^3 and /X4; 
however, we will keep them sufficiently small (say O(lGeV)) that we only need consider 
perturbations about the positions of the minima computed in section 3. We must return 
to ( |3.4|) and calculate the corrected positions of the minima. Let's see what happens to the 
minimum in the (fix direction. It will pick up small components in the <p2 and ^2 directions. 
The first two equations in ( |3.4| ) may be combined for non zero <p\ to yield: 

2/^50102^2 - nl<h + (4to = • 

The first term here is second order in the small fields and so we obtain for the tangent 
of the CP violating phase: 1^2/ <l>2 ~ 1^1/ 1^3- Similarly, the other minimum will pick up 
a small component in the <j>\ direction. Either of the first two equations in fl3.4|) can be 
used to determine the magnitude of $2 by ignoring the small <fi\ terms. With non-zero <fii, 
however, we are no longer free to choose ^2 to vanish, and hence the CP violating phase 
to be zero, and instead must determine it from equation ( f4.1|) . Once again, the first term 
is negligible because <fi\ is perturbatively small and we obtain as before: 1^2/ <p2 ~ ^4/^3- 
The small component is then given approximately by the solution of the cubic derived 
from the third equation in ( |3.4| ), where, for <f>2 and 1^2 we use the values we have just 
calculated from the first pair of equations in (|3.4|) . 



The phase structure of a bubble wall can thus be estimated when we recall from the 
previous section that in the middle of the wall, where both </>i and </>2 are large, the phase 
tends to be driven to zero so as to minimise the k§ term in the potential. We can thus 
typically expect the variations in phase experienced by top quarks during reflection from 
the wall to be A9 = C^arctan/Ltl//^). Naively, in the spirit of 0, we might expect the 
preferential reflection of hypercharge to be proportional to A9 1 , which may lead, ultimately, 



1 We should note that, in the limit of the CP violating parameters tending to zero, the 
preferential reflection of top quarks with definite hypercharge will not be proportional to A9. This 
is a reason for wishing to obtain a description of the transition without making these parameters 
too small. 
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to a baryon asymmetry. Thus we see that the introduction of small 113 and /X4 values may 
result in net baryon production at the phase transition wall. 

For somewhat larger values of ^3 and ^4, we do not expect the essential features of 
the transition to be altered, but to obtain accurate results for the quantities involved, it 
becomes necessary to find the relevant solutions of (|3.4| ) by numerical methods. 

Recall from section 3, that at the jumping temperature, we have the following ratios 
of VEV to temperature : v/Tt = 0.67 in the false vacuum; v/Tt = 2.0 in the true vacuum. 
Baryon violation is suppressed in its rate by a factor of exp(— ^v/T) |I(J, where gw 
is the SU(2) L gauge coupling. When the VEV makes its jump during the second stage 
of the transition, the rate is thus reduced by a factor of approximately 4 x 10 -12 , to a 
value which is slow enough to ensure that any baryon asymmetry created will not be 
subsequently dissipated. In false vacuum regions, in contrast, the baryon violation, whilst 
not completely unsuppressed, can still be rapid enough to convert hypercharge to baryon 
number before the bubble wall sweeps past. To generate the observed baryon asymmetry, 
therefore, we favour bubble walls with low velocities and a large CP violating parameter 
A6 = (9(arctan/4/At|)- 

5. Conclusion 

In exploring the parameter space of two Higgs models, we have found that the 
symmetry breaking process may occur in two steps, the second of which produces many 
of the features of a strongly first order phase transition. We have shown that this may be 
important in certain weak scale baryogenesis scenarios, where it is crucial that the baryons 
do not disappear immediately subsequent to the transition, as happens if it is only weakly 
first order. 

Calculation of higher order thermal loop graphs to obtain an improved effective 
potential will only yield small corrections since the interesting stage of the transition 
occurs well away from the origin in field space. It may be interesting, however, to work 
out more accurately, just how large a CP violating phase can be obtained by adjusting the 
various parameters. 
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Figure Captions 



Fig. 1. Temperature development of the positions in field space of the minima of the 
effective potential and the vacuum expectation values of the scalars. 
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